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Every effort is made "by the NACA to insure that its published 
translations are accurate reproductions of the original work of the 
authors. Papers are selected for translation on the basis of interest 
and probable usefulness and, although an examination is made for 
technical soundness, the Committee cannot assume responsibility for 
the accuracy of the detailed results presented by the author in the 
original paper. The Committee will of course call attention to any 
errors observed at the time of publication or subsequent thereto. 

Analysis of NACA Technical Memorandum No. 1150 has disclosed 
that the principal result of the original paper is in error. The 
compatibility of the three branches of a X - shock (Y- shock, branch- 
shock, fork-shock) in a supersonic flow is granted at crossing points 
of two supersonic arcs corresponding to two opposite families of 
characteristics and one common subsonic arc ending at the sonic line. 

The further continuation is sometimes rather the characteristic of 
the same family than the supersonic arc of the shock polar itself 
corresponding to the opposite family. Such cases with characteristics 
are not considered by Eggink. Therefore, his classification of ^.-shocks 
and the sensational result that X - shocks cannot exist below Mach 
number 1.245 without any value. Obviously k-shocks are possible at 
all supersonic velocities and their appearance is without any special 
physical significance. 
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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL MEMORANDUM NO . 1150 


COMPRESSION SHOCKS OF DE 'ACHED FLOW* 

By Eggink 

SUMMARY 

It is known that compression shocks which lead from supersonic to 
subsonic ■velocity cause the flow to separate on impact on a rigid 
wall. Such shocks appear at. bodies with circular symmetry or wing 
profiles on locally exceeding sonic velocity, and in Laval nozzles 
with too high a back pressure. The form of the compression shocks 
observed therein is investigated. 


I . INTRODUCTION 


The present investigation of the forms of shock proceeded from 
schlieren observations of wing profiles and an investigation of the 
forms of shock in Laval nozzles with too high a back— pressure . 
Schlieren photographs of the shock forms in Laval nozzles have 
been published already in FB 1756 . 

In observing shock it is necessary to be sure to have a 
suitable adjustment of the schlieren optical system. Good 
definition of the shock form can be obtained with the schlieren 
edge perpendicular to the flow, and bright shock lines against 
a dark field. Small tunnel breadth yields sharp shock lines, 
while the shock is no longer perfectly plane and appears streaked. 
The observation is often disturbed by violent fluttering of the 
shock. Then the shock appears to the eye as a broad streaked band 
and in photographs, as a set of shocklines. (Compare fig. 3*) 


*"lJber Verdichtungsstbsse bei abgeloster StrUmung , " Zentrale 
fur wissenschaftliches Berichtswesen der Luftfahrtf orschung des 
Generalluft zeugma isters (ZWB) Berlin-Adlershof , Forschungsbericht 
Nr. I 85 O, Aachen, Aug. 12, 19^3 • 
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II. THE ASSORTED SHOOK FORMS 


It will be assumed that there is a compression shock in a 
supersonic floxv. If this shock encounters a rigid wall, there 
are three oases to distinguish. (See figs. 1(a) to 1(c).) 

(a) The shock runs perpendicular to the wall; the flow does 
not separate. 

( b) The shock runs oblique to the wall; the flow separates 
from the wall. 


(c) The shock branches into two separated shocks in front 
of the wall. 

In case (b) the flow generally returns to the wall after flowing 
past the dead space,- since the pressure gradient directed against 
the wall downstream from the bent shock forces the streamlines 
against the wall. From the observation that the flew behind the 
shock is along the wall again, it cannot be inferred, therefore, 
whether case (a) or (b) is present. (Compare fig. 2.) 

In contrast to case (b), which is observed shortly after 
exceeding sonic velocity, the branch of (c) (fig. 3) is first 
observed at higher Mach numbers. In agreement with this, it follows 
from the theory of the branching in Section III, that a branching 
first oomes' into existence for 


w 

M*= ~>1.191 

a* 

or 


M 


_ W; 

a 


1.2JU-5 


The schlieren observations of the branch show, that the 
principal shock A-D is bent, in general, and the branch shock 
lines A-B and A-C can be regarded as practically straight. By 
this the stream boundary B-C is taken as straight and the pressure 
is assumed constant in the dead space along this boundary. 

(If the flow is already detached weakly upstream by considerable, 
roughness of the wall, then a fan of shock lines, which come 
together in the principal shock forms in place of the branch.) 
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The flow in the vicinity of the branch point is commuted in 
Section III. This gives as a result the shocks coining together 
at the branch point A and, with that, the shock lines A-B and A-C 
as functions of the flow and the diversion of the flow in the 
separation point B or the pressure in the lead space. 


III. THEORY OF BRANCHING 
1. Statement and Execution of the Solution 


The individual compression shocks coming together in the 
point A are shown in the familiar shock-polar diagram of 
A. Busemann. The streamline A-E, which passes through the 
branch point, separates the two flows ( 3 ) and (J 4 .) (fig. !+) 
behind the shock. 

The pressure P in ( 3 ) and (1+) must be equal along this 
boundary, end the directions of the velocities in (3) and (4) 
be parallel to it. The different increase of entropy from 
( 1) to (1+) and (1) to (3) yields a discontinuity in the velocity, 
temperature, and density across the boundary. 

The flow ( 1) corresponds to the shock pole in the polar 

diagram (fig. I 4 .) on which ( 2 ) and ( I 4 .) lie. ( 3 ) lies on the shock 
pole Mg* which was defined by means of (2), which is variable. 

The conditions required at the separation line: 


p 3 = p ii 


£ 2 + ^3 


( 1 ) 


suffice to ascertain the points ( 3 ) and (J 4 .) for each Mg* or &g. 

In the dimensionless formulas for tho compression shock the 

to' the static pressure of the flow, 

. The crocked shock A-B from ( 1) to (2) 

he fechn. Berichten Bd,. 10 ( 191+3 ) . 

Heft 2, Prof. A. Weise DVL likewise w r orked with branch shocks 
and computed the special case of branching with a straight 
principal shock. 


pressure p is referred 

JL pressure ratio | 



“ ^ As can be seen in • 


Po 


b 
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reduces the static pressure to p 0 * On the shock polar M 2 * the 
pressure is therefore to be referred to ^ 0 . The abbreviation x 
is used for the pressure ratio of the flow that characterizes the 
shock polar, the abbreviation, y for the pressure ratio along 
the shock polar, and indicate with the pressure ratio referred 
to p c . According to figure i+s 


Elr-- 

Po 


XI 



Po 


P2 


= Y2 


P 3 

Po 


Ek 

Po 


= yk 


Equation (l) is referred bo p Q yields: 





y k 


b 2 + *3 ~ H (2) 

Po 

The ratio of the static pressures is known to be a function 

Po 

of y2 and x. The formulas for the compression shocks are 
assumed as known. 2 

The angle of deflection & of the flow follows as a function 
of y and :x from* 


2 For further detail concerning the fundamental equations of 
the compression shock and the formulas derived therefrom, see: 

R. Sauer "The ore t, Einfuhrung in die Gasdynamik, " Springer 19i|3* 
The designations chosen there have been retained here. 
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4o 

Po 


r" - «P^i.y2)® 



(to 


if equation (3) is set up with the same indices as in equation (2): 


f ( I i>5 , 2) t f (*2'-y 3 ') 


and with the use of equation (Z4.): 


ya _ y 2 

1 (MO ^(*1^2) 


as well as y? » - — from equation (2). 

fj>(xi,y 2 ) 


3 

Eliminating the shock angle a from equations (129) and (I30) 
on page 65 of "Theoret. Einfuhrung i.d. Gasdynamik," leads to 
equation (3)5 correspondingly, equation (k) is obtained' from 
equation (129) page 65 and equation (150) page 77. 
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y 2 ^4 

f(xi,yi^) = f(x lj ,y 2 ) + f — 

cp(^i^y2^ <P(xi,y 2 ) 


it is seen that the unknown y^ cannot he computed explicitly 
as a result of the complicated function f . 


The solution of equation (2) is obtained graphically. The 
following construction is suitable to this end (fig. 5). The 
distribution of the shock polar M x * (= MLj* - distribution) 

■8 = f(x 1 ,y) from equation (3) is plotted against the pressure 
ratio y. The d distribution of the shock polar Mg* 


■^o 

( = M-,* — distribution) is plotted against y — displaced 
^ Po 

parallelwise by The exaggeration of the y-axis by the 


8 


Jr O 

factor — refers the pressure of the shock polar M^* to p Q . 
P° 

With that, the starting point Xp’ , since x* — ** y 2 at the 

£ ^ P 0 


point * 2 , coincides with y g of the M3* - distribution. The 

points of intersection of the two distributions furnish the desired 
solution, since the conditions of equations (2) are satisfied here 

$ 

with $2 + = ^ and y^* p~ ” y^* 


In carrying out the construction, the $ distributions were 
computed for M x * = 1.1; 1.2; ... 2.2; 2.3 and for M 2 * at any 

instant the smaller value was chosen, for example, for Mj.* = 2.0, 

= 1,9; 1.8; . . . For M-.* there are intervals in which the 

C. ■*- 

solution shows various relations. The boundary values for these 
intervals as well as the "limiting points" mentioned in the 
following are treated in the next section. 


The manner in which the various solutions in the individual 
intervals are accounted for in terms of relation of intersection 
points is schematically illustrated in Section 6-8. 


If 1,191 > M3* > 1 then the M g * - distribution lies in the 
interior of the M-^* — distribution (fig. 6). A branching of the 
shock does not take place in this region, therefore. 

The interval 1.842 > M]_* > 1.191 . provides a point of intersection 
(fig. 7), At the starting point x^ the Mg* — distribution 
coincides with the M-|_* - distribution. The Mg* - distribution 
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displaced from x^ with respect to the distribution by a 

differential shows the "limiting point" a. The point of inter- 
section moves out from a with diminishing Mg* along the 

Mi* — distribution up to the limiting point b. On the variable 
M 2 * — distributions,, the intersection point approaches the starting 
point y 2 of the Mg* — distribution, in doing so. At the limiting 

point b, the intersection point coincides with the starting. point 
of the Mg* — distribution. Here Mg* reaches a lower limit below 

which there is no solution, since the Mg* — distribution lies within 

the Mi* — distribution for smaller values of Mg*. 

If Mi* > 1.842, there are two points of intersection, which 

run out from the limiting points a and c (fig. 8). At the 
limiting point a, the variation of the one intersection point 
known from the preceding interval commences, and terminates at the 
limiting point b. The second intersection point starts out at 
the limiting point c and coincides at the limiting point d 
with the starting point of the appropriate Mg* — distribution and, 

thereby, terminates the ambiguity. Near a lower limit for Mg*, 

below which there is no longer a solution (limiting point b) 
there is a lower limit at which the ambiguity vanishes (limiting 
point d) . 


2. Obtaining the Limiting Values Limiting Points a and c 

At the intersection point of the M-j* - distribution and A 

. ; V p o 

Mp* — distribution, displaced by Ax, d» = f(x^ + Ax) y T ~~ 
^ Po Po 

the following is valid' according to figure 9s 

a a d* + y* =- " y * 

Po g 


or 


and 


lim 
Ax— >0 


/\ 

, ' X p 0 

f Kxi + Ax) -- y 


A - 

P 

f(x x + Ax) — yg 
Po 


Ax 


- f(xi,y g ) = -Ad 


f(*l,y J 

x 8 - lim ^ 

- Ax 

Ax— ^0 


I 


I 


mu 


mi 


■ III 


3 BACA rn Bo. 1150 

and right side goes to tan £i in the transition at the limit, 
according to figure 9 5 this is the direction of the Mi* - 
distribution a.t xj. • 

It is: 


tan pi = 


df (xj,y) 

dy 

y 


*1 


the left side goes over the nartial derivative of f with respect 

P' 

to x at x = xi since --- ■ —»> '! for Ax — > 0 . It follows that 

p o 

(x,yg) 

5 x 

X = Xl 



The quantity y g which defines the limiting points a and c is 
computed from equation ( 5 ). Differentiation of equation (3) yields: 



6 : 


y = x. 


t 


2 - (k + 1) xi 


k-1 

IT 


k + 1 

xi 2k 


a - ki 


k-l 

' k 


tan pi 


2k 


( 6 ) 
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in which $ is to be taken from equation (3) with x = and 

y ~ y g ,* 

The quantity y g, defining the limiting points a and o 
can be obtained only by aporoximation from equations (5)» (6), 
and ( 7 ) • 

Limiting Points b and d 

According to figures (7) and (8) the .'intersection point of 
the Mq* - and the M 2 * - distribution falls on the limiting 


10 
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points b or d at the starting point of the Mg* - distribution. 
Here, therefore, both curves have two points in common, that is, 
they have a common tangent, or the direction of the M]_ - • 

distribution = tan a.2 at the point yg agrees with the direction 
of the Mg* - distribution at ' the starting point = tan Pg. 


tan 02 = tan p 2 ( 


The direction at the starting point of the Mg* - distribution 
follows from differentiation of equation (3)* (Compare equation (6).) 



In the denominator the exaggeration of the 


factor 
be taken 


, mentioned above' 
Fo 

into consideration. 


for the Mg* - 
If 


y-axis by the 
distribution, must 




P2 1 1 

P °'. to/Po ^ ?Y P ° 



and 


Po 


from equation (1+) are substituted, it follows after calculation*. 


tan P2 = 



2kj/k - 1 


(k - 




x 


1 



k - 1 
k + 1 




( 9 ) 


Differentiation of equation (3) gives the direction of the M]*/- distribution at the 
point y = y 2 with 3 from equation (3) as: 


tan a 2 = 


1 sin 2 -d 

2 



1 


k -fr 1 
k - 1 


X 1 





1 



k - 1 
k + 1 


X 1 



( 10 ) 
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Here too, the quantity yg sought for is obtained only as an 
approximation from equations (8), ( 9 ), and (10). 


The Interval Limits 

Figures 10 to 12 show the variation of the curves tan 
and tan £2 with yq for fixed Hq*, in the various intervals. 

The intersection points of both curves furnish the yg - values 
of the limiting points b and d. 

Both intersection points approach the starting point xq of 
both curves with diminishing Ml*. First of all, d, coincides 
with xq j because of this ambiguity vanishes, since; only the 
limiting point b still exists. The variation in curve of 
figure 10 goes over to the variation of -the curve in figure 11. 

Then b coincides with xq and the variation of the curve 
corresponds to figure 12. The vanishing of the limiting point b 
is only possible, however, if the M 2 * - distribution lies within 
the Mp* - distribution. Here the lower limit for Mq* is reached, 
accordingly, below which there is no solution. The .curves tan 
and tan £2 agree at the point y = xq in magnitude and direction, 
that is: 

tan ~ "ten P2 (H) 

b tan cc 2 _ 6 tan ^ ^ 2 ) 

by by 

The somewhat detailed proof of these relations by equating and 
differentiating equations (9) and (10) is passed over here. 

Since the directions of both curves at the point Xq agree, 
the coincidence of the intersection points b and d with xq 
means that the curvature. of both curves agree at the point xq. 
Therefore, with equation (12) the two second derivatives of both 
curves are equal. 
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d 2 tan 

b y 2 



tan Pg 



(13) 


The interval limits sought for follow from, equation ( 13 ) . 

Taking the partial derivative of equations (9) and ( 10) 
and substituting in equation (lj) leads after a detailed calculation, 
which will not he entered into here, to an algebraic equation of 
the fourth degree in x^* By means of the solutions x^ a =1 

and x-j^ ss 1, these reduce to the following, equation of the 

second degree: ‘ ■ 




3k - 3 
k + 1 


0 


with the solution: 


k-1 

k 


(k - 2) ± ,/ k-L-i. 

'J k ♦ 1 


or with: 


xi 


k-1 

k 


1 - 


k-1 
k + 1 


M 


r *2 


-1 


solved for 


Mj* 2 = (k + 1) ( 1 + 


{ \±zJL 

If k + 1 


This gives: The lower limit for the ambiguity lies at - 1.81+2 

or ¥]_ ~ 2.569? the limit, below which there is no solution, lies 


at M* = 1.191 or M x — 1.21+5. 
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3 . Plotting the Polar-Shock Diagram 

\ 

The result of the construction and calculation is displayed 
in the polar-shock diagram (see appendix). 

is given, with calculated value for y l±, according to 
figure^ and the . point (3) ■ according to figure If,. The points ( 3 ) 9 -re 
computed for the flow' at any given instant are plotted 

on the variable -shock polar s Mg*, and joined (fig. 13) by a 
curve. The curves begin at the shock polar Mq*, at the limiting 
point a or c and end on the axi s of the diagram at M* that 
corresponds to the limiting points b and d. In the shock- 
polar diagram the point (Ij.) or & ^ (fig. I 4 .) , for given point (3) 
or ■d- , is defined from equation (’ 2 ) t 

^ *' -t . 

^ Lf~ ^ 2 + ^3 


7/hether (U) lies on the supersonic or the subsonic side of the 
shock-polar follows from the variation of additional points. The 
lower limiting line for Mo*, below which there is no branching, 
is the connecting line of the limiting- points b (curve b). 

The lower limiting line for the ambiguity is the connecting of 
the limiting points d (curve d). 

The values of Mg* for which ^ = 0, therefore, for which 
the shock A-D becomes straight, are determined from the inter- 
section point of the curve e marked with little lines with the 
shock polar M^*. If e coincides with the starting point of the 
shock polar, then the limiting point a must lie on thq axis of 
tho diagram, since the limiting point a corresponds to the 
"solution" at the starting point of the shock-polar and 0, 

This condition used in equation ( 5 ), (6), and (7) yields: 



as the intersection point of the curve e with the axis of the 
'diagram. (It is know that the maximum of the pressure for this 
value is behind 9 . straight shock.) 

In figure II 4 . the branching with straight principal shock is 
shown for various , in figure 15 the branching with fixed 
Mq* = 1.5 and $2 variable. 
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IV. SUMMARY 


Schlieren observations show a branching of the compression 
shock in the vicinity of a fixed wall. The flow in the vicinity 
of the branching point is computed, and the result displayed in 
a shock-polar diagram. 

A branching of the shock only appears for Ml* > 1.191 
(Mi > 1,21+5)* 2? Mi* > 1.81+2 {Ml ^ 2.5U6), then two solutions 

exist. If Mi* lies between 1.191 and 1.355 the principal 
shock is oblique, the branching deflects the flow direction 
positively, only for Mi* > 1.353 (Ml > 1.1+81+) the flow direction 
is deflected both positively and negatively, here a branching with 
a straight principal shock becomes possible. The second solution 
for Mi* > 1.81+2 yields only a positive deflection of the 
flow direction. 


Translated by Dave Feingold 
National Advisory 
Coxnmitt"!© for Aeronautics 
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Fig. 9 



9. Determination of the limiting points a and c. 
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Figs. 10-12 




Figures 10-12 Determination of interval limits. 
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